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Abstract 

We consider directed percolation processes for particle types A and B coupled 
unidirectionally by a transmutation reaction A — > B. It is shown that the strong 
coupling regime of this recently introduced problem defines a universality class with 
upper critical dimension d c = 6. Exact expressions are derived for the scaling di- 
mensions in the inactive phase above d = 4. Below d = 4 the interactions of the 
normal directed percolation also get relevant. 

PACS numbers: 64.60.Ak, 64.60.Ht,64.60.Kw. 

It is a well known fact that the phase diagram (or more generally, the parameter space) 
of homogeneous media may contain critical submanifolds. Near such a critical manifold 
some local restoring forces vanish and the system becomes "soft" with respect to some 
external forces. As a consequence next-neighbour interactions and long-range correlations 
begin to play an important role. If in addition a random force is acting in the system, 
then a complex behaviour may result. The most important experimental aspects of critical 
phenomena are scale invariance and universality. 

The universality classes may be subdivided into categories like "Critical Statics" , "Crit- 
ical Dynamics" , "Reaction-Diffusion" and "Self-Organized" . Dozens of universality classes 
have been identified. Formally there are infinitely many universality classes - most of them 
without a plausible physical interpretation, of course. 

The directed percolation (DP) universality classpQ[2] is one of the simplest and most 
important in the reaction-diffusion type category. It occurs in contexts like birth-death 
processes, population growth and chemical reactions. [3] It also may be given a purely 
static interpretation |Tj[2] and it formally coincides with Reggeon field theory, jl] This field 
theory contains nonlinearities of lowest (=cubic) order, and the general expectation is 
that it thus describes various processes with an active state for a one component order 
parameter. Essential preconditions are that the particles diffuse and that there are no 
conservation laws. 

Recently Tauber et al. [S] introduced a model describing DP processes for N particle 
types A, B, ... coupled by transmutation reactions. For 2 particle species the model may 
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be defined by the reaction equations 
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2 A, A ► 0, 2 A ► 0, A n .i, 1+e . 
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B -> 2B, B^O, 2B^0, B n A B n+e , 
A^B. 

The particles may branch, decay and annihilate with some given rates. Together with the 
diffusion reaction (rate A) this yields normal directed percolation. The two DP processes 
are coupled by the transmutation reaction (rate a). An interesting interpretation of the 
reactions (0) arises when one considers the particle species as particles in different layers 
on a substrate. The reactions (JTJ) then are a model for the adsorption and desorption of 
particles. jH] 

A field theory corresponding to the process (JTJ) may be derived by standard techniquesjn] 
[7][H], the action functional reads [5 J 
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Here <fi corresponds to A, if) corresponds to B, and quartic terms that are irrelevant in the 
renormalization group sense have been dropped. The constant is 2 l ~ d T[~ d l 2 /V (d/2). To 
reach the critical point the parameters 7q and Tq (or their renormalized counterparts) 
have to be adjusted to 0. Here we shall only consider the cases r$ 
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For a transmutation rate a = S decomposes into uncoupled A and B DP action 
functionals. In fact, it follows from the reaction equations or the Langevin equations 

corresponding to S, that response functions of the type ^0 m 0"^ and (^p m ip n ^ are a- 

independent. For A the reaction A — > B has the same effect as A — > 0, and for B 
a has no effect as long as there there are no A particles. The coupling is only visible 

in the response functions {^) m <^j. This implies that the critical exponents /3, v and z 
retain their DP values. [H] It is also of interest to note that g^ and ux may be given the 
same value by rescaling <fr or <fr, making the action invariant under the transformation 
t < — > —t [ in the a = case. It is equally well possible to fix gx = 1 



and to get a renormalization group fixed point = 0(1) and — O (4 — d) instead of 
#0 = a/O (4 — d) and = \J O (4 — d), and analogeously for g^ and u$. 

Tauber et al. [5] examined the model (J2J) for a ^ and noticed that additional cubic 

interactions f d d xdt jsi^V'^ + ^iptfi 2 + ifV'V + s^ipcj)^ are generated in the loop expan- 
sion and must be taken into account in a renormalization group calculation. The complete 
model displays two fixed lines in the space of the coupling constants, one of them stable. 
This result was confirmed and generalized by Janssen||5] , who also showed that the param- 
eter a carries the same scaling dimension as tq — tq c . At any rate, a is strongly relevant, 
and a quick crossover to a strong coupling regime will occur. 

Here we are interested in this strong coupling regime and now set a — 1. For d / 
there is a "rapidity" symmetry [S] 
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Dimension 



> 12 

> 6 

> 8 

> 6 

> 8 

> 8 



Table 1: Field exponents and dimensions of possibly relevant terms (for the Ui = model). 
The asterisk denotes an exponent > 1. 



In a formal sense the term unidirectional coupling thus is somwhat misleading. To simplify 
the calculation we rescale the fields to get g = = g^. 

Dimensional analysis 

The goal now is a 2-loop renormalization group (RG) calculation for the model (J2J) in the 
<t = 1 case. A first and elementary step of any RG calculation is a dimensional analysis 
of the action integral. From this one obtains the critical dimension of the model and the 
canonical wave vector dimensions of the fields, coordinates and coupling constants. 

There are some choices here, and it is of interest to formalize the (simple) procedure. 
The condition that the terms of the action integral be dimensionless leads to linear equa- 
tions for the canonical wave vector dimensions of the fields and coordinates. Let us call 
the number of fields plus the number of coordinates the model order of the field theory (6 
in the case of Eq.Q). The determination of the model order minus 1 canonical dimen- 
sions of the fields and coordinates and of one coupling constant (or the critical dimension 
instead) requires at least model order linearly independent terms from the action integral. 
For e = the linear equations may also be interpreted as equations for a hyperplane in 
a model order dimensional exponent space. The integer valued coordinates of this space 
are the exponents occurring in the terms of the action integral, the canonical wave vector 
dimensions are the components of the normal vector of this hyperplane. Selecting the 
terms of a critical model thus amounts to determining a hyperplane in a model order 
dimensional space. The normal vector of this hyperplane (1, 2; 2, 2, 4, 4) in the case of 
Eq.(J2J)) together with the critical dimension, is a good signature for a critical model. 1 The 
terms in the halvespace containing the origin are relevant, the terms in the other halfspace 
are irrelevant. The procedure of selecting a critical hyperplane may also be captured in a 
computer program. 

For the action integral (j2J) the result for a hyperplane containing the if)(p-term is a 
critical dimension d r = 6, and with e = d r — d, 



[9] = e/2, 



Ui 



-2 + 6/2, [« 



2 - e/2, 



[V]=4-e/2. 



Since it is crucial for the consistency of the model it must be checked now that the action 
(j2*j) with Ui = and without the Sj terms contains all terms relevant or marginal near 
d = 6, that may be generated from other terms present originally, even irrelevant ones. 
The essential fact to note to this purpose is that any vertex with an external <fi at least also 



lr The normal vector of the critical hyperplane for a = is (1, 2; 2, 2, 2, 2). 



3 




Figure 1: Contributions to the T a vertex up to 2 loop order. On the top right there are 
also shown the graphs for the interactions. The heavy line symbolizes the if), the normal 
line the field. The dot denotes the r CT vertex in tree approximation. 

contains one external 0, and vice versa for if>: an external if? implies at least one external 
if?. This comes about because in the direction of increasing time a line may be converted 
to 0, but it cannot disappear. And in the direction of decreasing time a if? line may be 
converted to if?, but cannot disappear (see also Fig.(£Q)). In Table ((TJ we have listed the 
and if? exponents of potentially dangerous terms for combinations of and if? exponents. 
The consistency of the action (j2J) follows from the fact that the canonical dimensions of 
the field monomials are > 8 = d c + [u] , with the exception of the r terms that must be 
adjusted to anyway. 



RG calculation 

It follows from the dimensional analysis that the Ui couplings are strongly irrelevant and 
may be dropped, at least close to d — 6. The Feynman graphs for a 2- loop RG calculation 
for the simplified action (j2J) are depicted in Fig.Q. It is only the r CT vertex that gets 
renormalized. To generate a contribution to the interaction or to the propagator 
would require an irrelevant (fxf> m ifj n vertex in the left corner of the time ordered diagram, 
and analogeously for the interaction (rapidity inversion symmetry). Consequently only 

/ \l/2 

one independent Z factor Z a = ( Z^Z-i j is required. For the individual fields this leads 
to 

^ = ^ = Z (ji z^ = Z^, = Z a . 
The renormalized vertex function r CT then is 

Yf ) {k\ u, r, g R , n) = Z a T a {k\ u, r, g) , (4) 

where r = 7g — = t r and gn = Z a l ^ 2 g^T e / 2 is the dimensionless renormalized 
coupling constant. The condition that gn reaches a fixed point asymtotically requires 
Z (j ~ 

The evaluation of the graphs of Fig.(£Q) with standard techniques p[| yields 

A = pL 

B = n 
C = p 
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Using an external frequency — iu = fi 2 and minimal subtraction we then find the stable 
RG fixed point 

i&MinSu, = &r + (4 + 40 In (4/3)) e 2 + O (e 3 ) . (5) 

The e-expansion (j5J) does not seem to converge for e — 1, and it seems to be impossible 
to extract a numerical value for gff for higher e values. However, the critical exponents 
may be given in closed form. From Eq. (0J) and the /i-independence of the unrenormalized 
vertex function T a it follows 

r?> (-, - 2 , % g R m) = (/0 (-, - 2) \ g R (»)) ~ r. (6) 

This finally leads to the scaling equivalence relations for the response function, 

($<f^ (k,u,r) = (ipip\ T CT (j)<P) ~ l/k 4 - £ ~ l/u 2 - £/2 ~ l/r 2 - £/2 . 

This exact result is valid fur d > 4 and in the inactive phase r > 0. 

It was argued above that the response functions (^(f> m (j) n ^ and (^tp m tp n ^ are independent 

of a, and one might think that it is inconsistent with this normal DP behaviour to have 
nontrivial Z factors and anomalous dimensions for 0, 0, ip and if). Actually this is not the 
case, because general response functions of this type cannot be generated from the action 
with Ui = 0. 



Summary 

It was shown that the problem of unidirectionally coupled directed percolation may be 
solved with RG techniques above d = 4 by expanding around d c = 6. In the inactive 
phase (r > 0) the scaling dimensions may be given exactly. 

In principle the solution g = 0(e), Ui = found for 4 < d < d c may extended 

contineously to d < 4 with g = O ((4 — d)°) and U{ = O (4 — d). In this way the ^0 m n ^ 
and (^tp m ip n ^ response functions can retain their normal DP behaviour. This also agrees 

with the fact that the scaling relations <fi ~ if} ~ k 2 and ~ if> ~ k 4 ~ £ coincide with the 
uncoupled DP scaling relations at d = 4. 

The situation below d = 4 is more complicated for the (ip(fij response function, because 

here also the U{ interactions and possibly the mixed interactions considered by Tauber et 
al. [Sj get relevant. Tauber et al. derive a critical exponent fa = 1/2 — (4 — d)/8 + (4 — d) 2 
describing the expectation value of <fi in the active phase in an expansion around d = 
4. This procedure was criticized by Janssenj^] in that it relies on the assymption that 
an exponentiation of logarithms were possible. At any rate, for a finite coupling a the 
diagrams of the type shown in Fig.(£Q) must be taken into account - all of them are IR 
singular already below d c = 6. The situation also gets more involved in the active phase, 
where the irrelevant coupling constant plays an important role, limiting the growth of 
if): is a "dangerous irrelevant parameter" above d — 4, that becomes a normal relevant 
parameter below d = 4. 

That the 1-loop diagram " A" of Fig.(£Q) poses a problem for an RG calculation and 
an expansion around d — 4 was already noticed by Goldschmidt [TT] (diagram "c" of this 
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letter) and Tauber et al. 5j. An inspection of the diagrams for the model with Ui = 
indeed shows that /?2 is non-classical 1/2) already below d — 6. 

There is some similarity with the critical dynamics of the Heisenberg ferromagnet 
(model .7)^2], where the dynamics is nonclassical below d — 6, while the statics is non- 
classical below d = 4. Likewise the static critical exponents of model J are independent of 
the dynamics. An essential difference is that the model J dynamics adds new coordinates 
and fields, leaving room for a new critical exponent. In contrast, the scaling dimensions of 
the coupled DP problem are completely determined by normal DP, and one would expect 

A dimensional analysis for the general coupled DP problem with n particle species (or 
n layers) shows that here also only two nonlinearities of the type displayed in Fig.(£Q) are 
relevant near the critical dimension, which now is d c (n) = 2 (n + 1). These nonlinearities 
act in the first and the last layer, and below d c (n) more and more other nonlinearities get 
relevant. This indicates that the RG in the form used here may not be the appropriate 
technique to solve the general case. As a first step it would be of interest to fully understand 
the n = 2 case in the active phase and below d = 4. 
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